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_, ■ Abstract 
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>: 

O ! We present a new prescription for the resummation of contributions due to soft gluon emission 
QQ I to the transverse momentum distribution of processes such as Drell-Yan production in hadronic 
cn ' coUisions. We show that famihar difficulties in obtaining resummed results as a function of 
t^^ ■ transverse momentum starting from impact-parameter space resummation are related to the di- 
^ . vergence of the perturbative expansion of the momentum-space result. We construct a resummed 
Q ! expression by Borel resummation of this divergent series, removing the divergence in the Borel 
j> I inversion through the inclusion of a suitable higher twist term. The ensuing resummation pre- 
scription is free of numerical instabilities, is stable upon the inclusion of subleading terms, and 
the original divergent perturbative series is asymptotic to it. We compare our results to those 
obtained using alternative prescriptions, and discuss the ambiguities related to the resummation 
procedure. 
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1 Transverse momentum resummation 

The computation of transverse momentum distributions of heavy systems (such as dileptons, 
vectors bosons, Higgs) plays an important role in collider phenomenology, from the Tevatron to 
the LHC [1,2]. As is well known, the perturbative QCD expansion of the inclusive distribution 
contains to all orders powers of aslyi^igr/Q), due to the emission of soft and collinear gluons. 
When the transverse momentum Qj. is much smaller than the mass of the final state Q these logs 
become large and must be resummed in order for perturbative predictions to remain reliable. 

The resummation, to given logarithmic accuracy, can be performed [3] for the Fourier trans- 
form of the differential cross-section -^ with respect to Qt- Upon Fourier transformation, g^ 

turns into its Fourier conjugate, the impact parameter b, and large logs of Qt/Q become large 
logs of bQ. Fourier transformation is necessary in order for the contributions included by resum- 
mation to respect transverse momentum conservation, thereby avoiding the spurious factorial 
growth of resummed coefficients [4]. However, the Fourier transform must be inverted in order to 
obtain resummed predictions for physical observables. This is problematic because the Fourier 
inversion integral necessarily involves an integration over the region of impact parameters where 
the strong coupling is not well defined because of the Landau pole. 

This problem has been treated with various prescriptions. One possibility is to modify 
the behaviour of the strong coupling in the infrared in the Fourier inversion integral [3] {b^, 
prescription, henceforth): this procedure is widely used, but it is known to lead to numerical 
instabilities when the resummed results are matched to fixed-order ones [5]. A second option is 
based on the observation that the Fourier inversion integral can be computed order by order in an 
expansion of the resummed results in powers of as'- if only leading log terms are retained in the 
Fourier inversion, the result is then well defined for all values of Qt [5]. This procedure however 
is unstable to the inclusion of subleading corrections: the Fourier inversion can be performed 
to next-to-leading log accuracy [6] (as it is necessary if the resummation is performed to this 
order), but in such case the result differs significantly from the leading log one, and in fact for Q 
around 100 GeV it blows up for values of Qt of order of several GeV, well within the perturbative 
region. A "minimal" prescription which is free of these difficulties can be constructed [7] , along 
the lines of the similar prescription for threshold resummation [4] . Namely, the integration path 
in the Fourier inversion is deformed in such a way as to leave unchanged the result to any finite 
perturbative order, but avoiding the Landau pole and associate cut in the resummed result. 
This leads to a prescription which is free of numerical and perturbative instabilities: its only 
shortcoming is that it is difficult to assess the ambiguities related to the resummation procedure, 
as it can be done in the 6^, prescription by varying the way in which the infrared behaviour of 
the strong coupling is modiffed. 

Here we shall show that, analogously to what happens in the case of threshold resumma- 
tion [8], the ambiguity in the resummation procedure is due to the fact that the perturbative 
expansion of the resummed result for the transverse momentum distribution itself in powers 
of as diverges. After discussing, in the next section of this paper, how existing prescriptions 
treat this divergence, we will show in section 3 that the divergent series can be treated by Borel 
summation, as is the case for threshold resummation [8,9]. The Borel transform of the series 
converges and can be summed. The inversion integral which gives back the original series di- 



verges, but the divergence can be removed by including a suitable higher twist term. This leads 
to a resummed result of which the original divergent series is an asymptotic expansion. The 
ensuing prescription is given in terms of a contour integral which is easily amenable to numerical 
implementation. The result is free of numerical instabilities, and stable upon the inclusion of 
subleading corrections. An estimate of the ambiguity on the resummed results may be obtained 
from a variation of the higher-twist term which is included in order to render the results conver- 
gent. In section 4 we will compare the result of our prescription to other existing prescriptions 
in the case of the Drell-Yan process, and discuss the ambiguities related to the resummation 
procedure. Some results on Fourier transforms are collected in the Appendix. 

2 The need for a resummation prescription 



Let us consider a parton-level quantity S which depends on a large scale Q and a transverse 
momentum g^,, such as the partonic Drell-Yan differential cross-section J^. Resummation is 

necessary because the perturbative coefficient of order n in the expansion of S in powers of 
as(Q^) has the form 
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P„(lngy) is a polynomial of degree 2n — 1 in Ing^, Qn{(jl) is regular as g^ -^ 0, and Dn are 
constants (see the Appendix for a definition of the + distribution). Physical observables are 
obtained, exploiting collinear factorization, as the convolution of parton level cross-sections with 
parton distributions [3] . When Q^ is large enough, it sets the scale of parton distributions, and 
the Qt dependence is entirely given by the partonic cross-section. For lower values of Q^ the 
scale of parton distributions is set by the impact parameter b, which is Fourier conjugate to q^, 
the convolution must be performed in b space, and the Fourier transform must be inverted to 
obtain physical predictions. In either case, the resummation is performed in b space at the level 
of partonic observables. 

Upon Fourier transformation, qj. is replaced by its Fourier-conjugate variable, the impact 
parameter b, and the small-g^ region is mapped onto the large-6 region. Large logs of b can then 
be resummed, leading to an expression of the form 
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is the large logarithm which is resummed, and 0{L^) denotes terms which are not logarithmically 
enhanced as 6 ^ oo. For future convenience, we have introduced in the definition of L an 
arbitrary constant bg (to be discussed below), and we have further defined 

a = /3oa.(Q'), (2.6) 

Po is the first coefficient of the QCD beta function, 

Q' ^^^ = -/?o aliQ') [1 + A asiQ') + OiaD] (2.7) 
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The inverse Fourier transform of S with respect to b is given by 
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using two-dimensional polar coordinates for b = bQ, and the integral representation of the 0-th 
order Bessel function, 

Jq{z) = — [ d9e-'''°''^. (2.10) 
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Now consider specifically the resummation of 

S(a„g2) = li^, (2.11) 

where ^ is the partonic transverse momentum distribution of a massive final state, and ctq the 
Born-level total cross-section. In this case, the 6-space resummed result has the form [3] 

J2{as,aL) = exp S{as,d L), (2-12) 
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A{as) = Aias + A2al + ...] B{as) = BiOs + . . . , (2.14) 

and the constants Ai,Bi can be determined order by order by matching to the fixed-order 
calculation. 

The integral in eq. (12.131) can be performed explicitly, and the result can then be expanded 
as 

oo 

^(a„aL) = ^a*-^/-,(aL), (2.15) 

i=0 



where inclusion of the first k orders in the sum corresponds to the next^-to-leading log (N'^LL) 
approximation. The LL and NLL functions /o, /i are explicitly given by 
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Note that with y = aL, using the leading log form of agi^Q"^) 
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It is apparent from eqs. (I2.16f2.17p that Tj{as-,aL) has a branch cut along the negative real 
axis in the complex plane of the variable y = aL: 



Re (y) < -1; Im (y) = 0. 



(2.19) 



This is due to the fact that the strong coupling blows up when its argument reaches the Landau 
pole, so that S{as,(yL) eq. (I2.13P is singular when b becomes large enough, i.e. when 
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At leading order, b^ = j^- It follows that the series for J2{as,aL) eq. (12. 4p has a finite radius 
of convergence, and the integrand in eq. (12.91) is not analytic in the whole integration range 
< 6 < +CXD, so the Fourier inversion integral is not well-defined without a prescription to treat 
the singularity. 

As mentioned in the introduction, various prescriptions of this kind have been proposed. 
Before discussing them, let us show that the reason why a prescription is needed is the divergence 
of the expansion in powers of ^^(Q^) of the resummed result obtained computing the inverse 
Fourier transform eq. (12.90 with S(as,aL) eq. (12.120 . To any finite perturbative order, the 
gr-space resummed result is found by expanding eq. (12.120 and inverting the Fourier transform 
order by order: 
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where we have replaced the argument q'^ of S by 
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When K ^> oo the series eq. f l2.2ip diverges. To see this, we compute the integrals in 
eq. (12.211) using eq. flA.l|) of the Appendix: 
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where the function Miji) is defined in eq. ( ]A.2p . we have assumed q^ ^ 0, so that distributions 
can be ignored, and the term with j = k, which leads to a vanishing contribution to T^xios, L), 
has been included in the sum over j eq. (I2.24p for later convenience. We now change the order 
of summation, and use the identity 

" d^e^r'^'-^^-', (2.25) 
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where the integration path H is any closed contour which encloses the origin .^ = 0. We obtain 

R.(a.,L) ^ 2±^±^K^^L^-' (2,26) 

Because of the singularity eq. (12.191) , the power series eq. (12. 4p has a finite radius of conver- 
gence equal to one 
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which immediately implies the vanishing of the radius of convergence of the sum over k in 
eq. (12:271) . 

The situation is thus similar to that which is encountered in threshold resummation [4,8,9]: 
the resummation is performed on quantities which are related by Mellin transformation to the 
physical ones, but the resummed results cannot be expressed as a Mellin transform of some 
function. Namely, their inverse Mellin transform does not exist, as a consequence of the fact 
that the inverse Mellin transform of their expansion in powers of asiQ"^) diverges. In the present 
case, the divergence of the perturbative expansion implies that the Fourier inversion integral is 
ill-defined; of course the problem disappears if one retains only a finite number of terms in the 
resummed expansion [10,11]. Various commonly used prescriptions replace the ill-defined integral 
with a well defined one, as we now review. In the next section, we construct a prescription which 
is instead based on the idea of replacing the divergent series with a convergent one through the 
Borel summation method. In the last section we will compare the various prescriptions and in 
particular the way they treat the divergence of the perturbative series. 



In the prescription of ref. [3], the variable b is replaced by a function b-^{b) which approaches 
a finite limit fonm < b^ as b —>■ oo, such as for example 
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In this way, the cut eq. fl2.19p is never reached. This procedure has some degree of arbitrariness 
in the choice of the function b^{b), which is interpreted as a parametrization of non-perturbative 
effects, whose size can be estimated by varying b^, for instance by changing the value of 6iim- 
The matching of this prescription to the fixed-order result is however numerically unstable, as 
pointed out in ref. [5]. 

A different possibility [5] is based on the observation that if only the leading log contribution 
(i.e. the terms with j = 0) are included in eq. (I2.24p . then the series converges, and its sum can 
in fact be computed in closed form, with the result (see eq. (JA.16P of the Appendix) 



SLL(a„ L) = 2^S {as, al) . (2.30) 

Equation (I2.13P implies that S{as,aL) depends on 6^ through 0^(1/6^). Therefore, using 
eq. (12.181) . the LL expression eq. (12.301) is seen to become a function of as{q'^). Therefore, 
the leading log truncation of the perturbative expansion in powers of ^^(Q^) eq. (12.271) has a 
finite radius of convergence, set by the Landau pole 

ql>Q^exp(-^j =A^, (2.31) 

where the last equality holds at leading order. 

The main defect of this result is that it is subject to large next-to-leading log corrections. In 
fact, the NLL Fourier inversion integral can also be computed in closed form [6]. The result (given 
in eq. ( 1A.17I) ) differs sizably from the LL result even for relatively large values of Qt (several GeV 
for Q = 100 GeV), as we shall see explicitly in Sect. 4 below. In fact, it turns out that the NLL 
correction diverges at a value of q^ which is an increasing function of the scale Q. This instability 
can be understood as a consequence of the fact that the truncation of the resummed result to 
finite logarithmic accuracy leads to an expansion in powers of ^^(q't) with coefficients depending 
on ln(gr/Q), where higher powers of asiQr) correspond to higher logarithmic orders. Such an 
expansion is necessarily poorly behaved at low Qt, all the more so when the scale ratio Qt/Q is 
large. Performing the Fourier inversion to leading or next-to-leading logarithmic accuracy thus 
removes the divergence of the series eq. (I2.2ip : this is analogous to what is found in the case 
of threshold resummation, where it can be shown [8] that the divergence of resummed results is 
removed if the Mellin inversion is performed to any finite logarithmic accuracy. However, the 
ensuing results are then perturbatively unstable. 

A yet different way of treating the divergence has been proposed more recently in ref. [7] , along 
the lines of the so-called Minimal Prescription of threshold resummation [4] . The basic idea here 
is that to any finite perturbative order, when the divergent series is replaced by a finite sum, one 
may choose the integration path in such a way that it avoids the singularities which appear at the 



resummed level. The result of the Fourier (or respectively Mellin) inversion is then unchanged 
to any finite perturbative order, but it becomes finite at the resummed level. It can be further 
shown [4] that the divergent perturbative expansion of the resummed expression is asymptotic 
to the result obtained in this way. This prescription is widely used [2]: whereas in the case of 
threshold resummation it leads to dependence of resummed physical results on a kinematically 
unaccessible region (albeit by power-suppressed terms), in the case of transverse momentum 
resummation its only shortcoming is speed limitation in its numerical implementation. 

3 The Borel prescription 

We now turn to the construction of a prescription which extends to transverse momentum 
resummation the Borel prescription proposed in refs. [8, 9] for the resummation of threshold 
logarithms. The basic idea is to tackle directly the divergence of the series (12. 24112. 2 7p by summing 
it through the Borel method. 

To do this, we take the Borel transform of eq. (12.271) with respect to a. This amounts to the 
replacement a^ -^ w^^^ /{k — 1)!, where w is the Borel variable conjugate to a. We obtain 

where in comparison to eq. (12.271) we have rescaled the integration variable ^ — ;> L^, and we have 
included all terms with 1 < /;; < j — 1, which vanish upon contour integration. 
Both sums in eq. (13. ip are convergent a,s K ^>- oo. Indeed, 
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the last condition being due to the simple pole of M(^) at ^ = 1. Thus, 

R{w, L) = lim Rk{w, L) = - / ^ e^« M(0 -^ S f a„ ^) , (3.4) 
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provided the contour H is chosen so that 

w<|^|<l. (3.5) 

Since M{^) has no singularities on the negative real axis, and T,{as,w/C,) has a branch cut on 
the real ^ axis between —w and 1, the integration contour can now be deformed so that R{w, L) 
is well defined for all positive values of w (see fig. [1]). 

The original function eq. (I2.24p is recovered by inverting the Borel transform: 
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Figure 1: The integration contour H in eq. (JST 



The inversion integral is divergent at w -^ oo. This is easily seen by inspection of fig. [T] as ly 
becomes large, the branch cut extends to the left, and the integration contour is pushed towards 
large negative values of ^, where M(^) oscillates with a factorially growing amplitude. 
We regulate the integral by cutting it off at w = C . We thus get 
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which is the Borel prescription for transverse momentum resummation. The result can be 
equivalently rewritten by doing a partial integration as 
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which may be more convenient for numerical implementations in that it depends directly on 
the physical observable S, rather than its derivative. Equation (13. 8p . and its equivalent form 
eq. (13. 7p . are the main result of this paper. It is interesting to observe that if we integrate by 
parts before cutting off the integral, then the surface term vanish. We then end up with the 
alternative resummation 
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As we shall see shortly, this is an equally valid prescription. 

In order to see that this is a valid resummation prescription, consider the truncation to order 



K of eq. fl3.6p , namely 
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where 
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is the truncated gamma function. The difference between the original Rk^ois^L) eq. fl2.23p and 
its Borel resummation R'^{as,L) is 
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i?^(as, L; C) is seen to be power-suppressed at large Q^ (higher twist): cutting off the w integra- 
tion at w = C is equivalent to the inclusion of a higher twist term, which cancels the divergence 
of the resummed expression. Specifically, i?^(as, Z; C) is a twist-t contribution with 

t = 2(l + C), (3.14) 

the choice C = 1 corresponds to the inclusion of a twist-four term. Moreover, it is apparent 
from eq. (13.121) that 

<(«„!; C) ~ e-§, (3.15) 



which vanishes faster than any power of as sls as —>■ 0. It follows that the original divergent 
Rk{c(s,L) is an asymptotic expansion of the Borel-resummed result R^{as,L) eqs. (I3.8|3.7p . 

Furthermore, the alternative prescription R^ (a^, L) eq. (13. 9p differs from R'~^[as, L) eq. (13. 8p 
by the first term in square brackets in (13. Sp . which is a finite higher- twist contribution. Hence, 
the two prescriptions correspond to two inequivalent but equally acceptable regularizations of 
the divergent sum which differ by finite terms, and are both asymptotic sums of the divergent 
series. 

The main features of the Borel prescription can be appreciated by considering as an explicit 
example of a resummed quantity T,{as, aL) = 7ll(«s, c^L), with 

7ll(«.,«L) = ^^-^^, (3.16) 



and S'll(«s,«-^) given by eqs. (12.1 3112. ISp evaluated at the leading log level (12.161) . namely 

7LL(as,«-^) = -^ ln(l + a-^). (3.17) 
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Substituting this form of S(a;s, aL) in eq. (13.71) . the associate g-r-space physical observable com- 
puted with the Borel prescription is found to be 

The C, integral is easy to calculate, because the integrand has only a simple pole at ,^ = —w: 

^^^{as,L) = —^- dwi^) M{-w), (3.19) 



where we have used the leading-log expression of the running coupling. It is thus clear that the 
divergent integration is cut off by the inclusion of a power-suppressed contribution 
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Note that the suppression is by powers of ^-: at finite order K the higher twist contribution is 
suppressed by a power of ^, as shown in eq. (I3.12p . but when resummed to all orders, the scale 



Q is replaced by an effective scale q: 
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4 Comparison of resummation prescriptions 

Let us now compare the results found using the Borel prescription to those of other prescriptions, 
with the dual goal of understanding the advantages and disadvantages of various methods, and 
of assessing the ambiguity which is intrinsic to the resummation of a divergent expansion. 

First, we look at a typical resummed observable. Namely, we consider the transverse mo- 
mentum distribution of Drell-Yan pairs, eq. (12. lip , which we evaluate at the partonic resummed 
next-to-leading log level, i.e. using eq. (12.120 with S{as,aL) computed including the first two 
terms in eq. (I2.15p . given in eqs. (I2.16|2.17j) with [12, 13] 
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Figure 2: The NLL partonic resuinined Drell-Yan transverse momentum distribution computed 
with various resummation prescriptions with Q^ = 10^ GeV^ and in a narrow (left) and wide 
(right) range of q^. 

The results are displayed in fig. [21 for Q = 100 GeV. The two lower curves at large q^ in this 
figure correspond to those found using respectively eqs. (1A.16P and eqs. (JA.17P of the Appendix, 
namely, to inverting the Fourier transform to leading and next-to-leading log accuracy (with 
bo = 2e'^""^). The sizable difference between these two results even for gj, as large as 10 GeV 
shows the instability of the truncation of the Fourier transform to finite log accuracy discussed 
in the introduction and first stressed in ref. [6]. 

The other prescriptions displayed in fig. [2] are the hi, prescription, where the Fourier inversion 
is performed after replacing h with 6^ eq. (12.291) . with h\\^ = bi, where bi = 7.2 GeV~^ is the NLO 
Landau pole eq. (I2.20p : the minimal prescription (MP) where the Fourier inversion is performed 
along the deformed path of ref. [7], and the Borel prescription eq. (13.71) with C = 1. 

In fig. [3] we further show the dependence of the Borel prescription on the parameter C 
which characterizes the higher twist term included in the resummation eqs. (I3.13f3.14p . as it is 
varied between twist four and twist eight. Because all these choices provide valid resummation 
prescriptions, this variation provides an estimate of the ambiguity which is intrinsic of the 
resummation procedure: indeed, the b^ and minimal prescription, also shown in this figure, 
are well within the band of variation as g^ ^ 0. These plots show that the ambiguity in the 
resummation procedure is negligible for gr ^ 5 GeV, it remains small for gr ^ 2 GeV, and it 
only blows up as g^ approaches the Landau pole. 

We can further elucidate the origin of these results by studying the effect of the various 
prescriptions when the divergent sum eq. (I2.2ip is truncated, so the Fourier inversion can be 
performed term by term. Consider specifically the first term in the series, namely, the inverse 
Fourier transform of L. The exact result is given by eq. flA.ip for k = 1, 
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The MP reproduces this exact result, because InlbQ/b"^) is analytic on the positive real b axis. 
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Figure 3: Dependence of the results shown in fig.[2]for the Borel prescription on the parameter C. 
The vertical line at q^ = 156 MeV indicates the position of the Landau pole. 



and a deformation of the integration contour has no effect; a branch cut on the positive real b 
axis only arises after summation of the whole series. 
The Borel prescription yields instead 
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as one can see by setting /ii = 4 and h^ = for all A; 7^ 1 in eq. (13.101) . The exact result 
is modified by the introduction of a correction of twist 2(1 + C). Note that the higher twist 
correction is tiny at large Q^, of order 10~^ for C = 1 and Q^ = 10^ GeV^. If we use the 
alternative Borel prescription i?*" (os,-^) eq. (13.91) we get instead 
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so the two prescriptions are indeed seen to differ by a higher twist term. 

Finally, the result of the replacement of b by 6^ eq. (12.291) can be computed analytically in 
terms of the Bessel function Ki : 



1 


1 bl 


1 
27r, 

2 
ql 



d^be-'^^-^\n 






62 



— I - biimqr Ki{biin,qj 



(4.7) 



Using the asymptotic behaviour Ki{z) ~ e ^ j \pz, we see that the correction factor in eq. (14. 7p 

z — *oo 

vanishes faster than any power of l/(6iimq'T) for g^ ^ l/^Um- 
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For higher order powers of L the same quahtative behaviour is found using the various 
prescriptions discussed here. Namely, the MP gives the exact Fourier transform eq. flA.ll) : the 
BP gives a result which differs from it by a higher twist term, and the fe^ prescription gives a 
result which differs from it by a term which is exponentially suppressed in l/(&iimq'T)- 

We thus see that the way different prescriptions tackle the divergence of the perturbative 
expansion is the following. In the LL and NLL case, the divergent series eq. f l2.24p is made 
convergent by truncating the Fourier inversion to finite order, i.e. by only retaining a finite 
number of terms in the inner sum over j. This, as discussed in Section 2, leads effectively to 
an expansion in powers of OLs{qT) which has very poor convergence properties at small q^ even 
when Q is large. The MP and BP both provide an asymptotic sum of the divergent series: 
the BP removes the divergence by inclusion of a higher twist term, and the MP by a suitable 
analytic continuation, which corresponds [4] to the inclusion of terms which are more suppressed 
than any power of Q^. At large Q^, the higher twist term of the BP is negligible so these two 
prescriptions are essentially indistinguishable when applied to convergent series. When applied 
to the divergent resummed expansion displayed in figs. [2}l3] they only differ in the region where Qt 
approaches the Landau pole, so the high-order behaviour of the series become relevant. Finally, 
the 6* prescription modifies the divergent series by inclusion of a term which is more suppressed 
than any power of l/(6iimgr). When applied to a convergent series, this prescription produces 
a result that differs sizably from that of the BP when q^ ^ Q^ and it approaches the Landau 
pole: this is because the scale of the correction term is set by Q^ for the BP, and by q^ for 
the h-i, prescription. At the resummed level, however, the effective scale of power suppressed 
terms becomes q^ also for the BP (compare eq. fl3.2UI) ). so all resummation prescriptions lead 
essentially to the same result. 

5 Summary 

We have constructed a resummation prescription for transverse momentum distributions which 
extends to this case the Borel prescription previously proposed for threshold resummation [8,9]. 
The construction is based on the observation that the reason why a resummation prescription 
is needed in the first place is that the perturbative expansion of resummed results in g^ space 
in powers of a^lQ^) diverges. The Borel prescription tackles this divergence by summing the 
convergent Borel transform of the divergent series, and then making the Borel inversion finite by 
inclusion of a higher twist term. The original divergent series is an asymptotic expansion of the 
result obtained thus. The Borel prescription is easily amenable to numerical implementation; 
being based on a 6-space resummation it is easy to match to fixed-order results, and it is 
perturbatively stable. 

There is some freedom in this prescription, parametrized by a real parameter C, related to 
the twist t of the term included in order to obtain convergence by t = 2(C + 1). Whereas C 
may be chosen to take any value, it is convenient to choose a value which corresponds to twists 
which already appear in the expansion of the observable being considered. Indeed, physical 
observables must be independent of the choice of C, and thus if an unphysical twist term is 
introduced, it must be compensated by an equal and opposite power suppressed term which is 
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thereby artificially introduced by this choice. 

Comparison of the Borel prescription to other available resummations, such as the minimal 
prescription or the b^ method, shows that at large Q^ they lead to results which are extremely 
stable and which only differ when q^ approaches the Landau pole. In fact, variation of the 
parameter C of the Borel prescription provides a reliable estimate of the ambiguity in the 
resummation procedure. For Q't ^ 2 GeV this ambiguity appears to be negligibly small, even in 
the region of a few GeV where the impact of the resummation is sizable. This is in contrast to 
the case of threshold resummation, where it was found [9] that the ambiguity is almost as large 
as the effect of the resummation itself in most of the kinematic region where the resummation 
is relevant. 

Our results contradict the widespread prejudice that transverse momentum resummation is 
affected by sizable ambiguities, and it shows that, at least as long as Q is as large as the W mass 
and Qt as large as the nucleon mass perturbative resummation of transverse momentum distri- 
butions provides reliable and stable results. The Borel prescription provides a new method for 
performing this resummation which has more stable matching properties than the 6* prescription 
and might be numerically advantageous over the widely used minimal prescription. 

Acknowledgements: We thank G. Altarelli for discussions. This work was partly sup- 
ported by the European network HEPTOOLS under contract MRTN-CT-2006-035505 and by a 
PRIN2006 grant (Italy). 



14 



A Appendix 

In this appendix, we collect some results on two-dimensional Fourier transforms of powers of 
logarithms. 

First, we compute the exact Fourier transform with respect to h of the fc-th power of In -^ 
(with 60 a constant). We get 



i- I d% e-'^^-^ In'^ ^ = 2 M(^) (0) 5(g^) + 2 5^ {^\ M^^^ ( 



0) 



— In'^-^' e 



where 

M(?7) 

and the + distributions are defined by 



hiyn^-ri) 



4/ r(l + r/)^ 



(A.l) 



(A.2) 



dql [D{ql)]^ = Q. 



In order to prove eq. ( JA.lj) . we define a generating function 



x(V")=(f)' 



V2 ak 



52 9?7* 



xikv) 



r)=0 



We have 



27r 



— d'be-'^-\{b,v) 



dbbJoibqr) (x^ 



(A.3) 



(A.4) 



(A.5) 



where we have used polar coordinates for b, and the integral representation of the 0-th order 
Bessel function 



Uz) = — dde 
The integral can be computed by means of the identity 



2tt 



-iz cos 9 



(A.6) 



+00 Y i- + - + -^ 1 

dxx^'JAax) = 2^'a-^-^^ 1 |f a > 0; -Rez/ - 1< Re/i < -. (A.7) 



We find 



We may now replace 



r(| + |-f) 



^ / d%e''^-''xib,v) = 2vMir^) [qiy' 



~2\^-l 



m' = m 



-2\^-l 






consistent with the definition eq. (1A.3I) . We get 
1 



27r 



d%e-'^^-\{b,ri) = 2M{r]) \ 5{ql) + 



_d_ 
dq'^ 



,2\^ 



m 



(A.8) 
(A.9) 

(A.IO) 
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Evaluating the k-th derivative of both sides with respect to 77 at r^ = leads immediately 
to eq. flA.l|) . Note that the term j = k is excluded from the sum because it vanishes upon 



differentiation with respect to q^. For q^ strictly larger than zero, both the term proportional 
to S{q'^) and the + prescription have no effect. 

Let us now turn to the evaluation of the Fourier transform to fixed logarithmic accuracy. 
Equation ( JA.ip shows that the Fourier transform of the k-th power of In 6 is proportional to 
l/g^ times the (fc — l)-th power of the log of the Fourier conjugate variable Ing^ (leading log 
approximation), but also includes terms proportional to all lower powers of this log. The N"LL 
approximation corresponds to including terms up to j = n in the sum in eq. (1A.ll) . i.e. such 
that the power of Ing^ is by ra + 1 units lower than the power of ln(6Q/6^). 

The NLL and N^LL approximations are particularly simple due to the fact that 

M«(0)=ln^ + 27s (A.ll) 

M(2)(0)= f\n^ + 2^E] (A. 12) 

where je ~ 0.5772 is the Euler constant. It follows in particular that if 60 = 2e~'^^, the NLL 
and NNLL terms in eq. flA.l|) vanish [12]. 

A useful form of the N"LL approximation can be obtained noting that 

M^^\0)= dxJi{x)W^. (A.13) 

Jo ^ 

It follows that eq. (lA.ip (for g^ > 0, i.e. neglecting distributions) can be rewritten as 



— I d%e-'^^-^\Yi^^ = 2—r t/xJi(x) Ing2+ln4 . (A.14) 

The N"LL approximation can then be obtained by retaining the first n terms in the binomial 
expansion of ( In g^ + In ^ j in this equation. 

This result is particularly useful in that it allows the computation in closed form of some 
Fourier transforms of generic functions to fixed logarithmic accuracy. Specifically, consider a 

function 

00 

F(L) = ^FfcL'=. (A.15) 

k=0 

Its Fourier transform to LL accuracy is given by 

1 /" - 1 d °° r°° d 

— d%e-'^^-'F{L) =2 — VFfc / dxJ^{x)\Yi^ql = '^T^F{\Tiql). (A.16) 

27rJ Jll ^^t^ Jo m 

This result was given in ref. [5]. 

One may think that because of eqs. flA.imA.T2|) eq. flA.161) with 60 = 2e~^^ automatically 



provides a result which is correct to N^LL accuracy. This, however, is not true if F{L) is 
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a physical observable, such as a cross-section. Indeed, in this case the N"LL approximation 
to it is defined by expansion of its logarithm: for example if F{L) is identified with J]{as,aL) 
eq. (I2.12p . the expansion of it to subsequent logarithmic order is given by the expansion eq. fl2.15p 
of S{as, aL) = In S(q;s, a L), and not of S(as, « L) itself. The NLL approximation to the Fourier 
inverse of F{L) may however be calculated exactly in terms of G{L) = In F{L). One finds 



1 
2^ 



d'be 



-iqj'-b 



F{L) 



NLL 



d 



dx Ji {x) exp 



Go + Gi In 



ft2 



d 



2 — F{\nql)M{G'{\nql)), 



(A.17) 



where 



with 



G{L) =\nF{L) = Go 



Go = G(lng, 



ft2 
Gi In ^ + O 

x^ 



Gi = G'(lng: 



In^^ 

x^ 



This is the result found in ref. [6]. 



(A.18) 



(A.19) 
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